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ABSTRACT 

Aims. The classic Weber-Davis model of the solar wind is reconsidered by incorporating alpha particles and by allowing the solar wind to flow 
out of the equatorial plane in an axisymmetrical configuration. 

Methods. In the ion momentum equations of the solar wind, the ion gyro-frequency is many orders of magnitude higher than any other 
frequency. This requires that the difference between proton and alpha velocity vectors be aligned with the background magnetic field. With the 
aid of this alignment condition, the governing equations of the multi-fluid solar wind are derived from the standard transport equations. The 
governing equations are numerically solved along a prescribed meridional magnetic field line located at colatitude 70° at lAU and a steady 
state fast solar wind solution is found. 

Results. A general analysis concludes, in agreement with the Weber-Davis model, that the magnetic field helps the coronal plasma to achieve 
an efi^ective corotation out to the Alfvenic radius, where the poloidal Alfvenic Mach number Mj equals unity (Mj is defined by equation <20» . 
The model computations show that, magnetic stresses predominate the angular momentum loss of the Sun. For the fast wind considered, the 
proton contribution to the angular momentum loss, which can be larger than the magnetic one, is almost completely canceled by the alpha 
particles that develop an azimuthal speed in the direction opposite to the solar rotation. The Poynting flux associated with the azimuthal 
components is negligible in the energy budget. However, the solar rotation can play some role in reducing the relative speed between alpha 
particles and protons for low latitude fast solar wind streams in interplanetary space. 
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1. Introduction 

The solar angular momentum loss rate £. consists of the par- 
ticle contribution and that contained in magnetic stresses 
Xm- The comparison of measurements of these quantities with 
models, the Weber-Davis analysis (1967 1 in particular, has 
yielded divergent results. Missions before Helios measured a 
total angular momentum flux X consistent with the Weber- 
Davis model (about lO-'^ dyne cm sr"'), but the measured az- 
imuthal angle of the bulk flow was generally greater than 1° 
at 1 AU (or equivalently 7 km s"' for an average slow wind 
of 400 km s ') (see Pizzo et al. 1983 and references therein). 
Such a large azimuthal flow speed implies that particles play a 
far more important role than magnetic stresses in reducing the 
angular momentum of the Sun. However, in the Weber-Davis 
model, 3/4 of the angular momentum flux at 1 AU is due to 
magnetic stresses. The Helios data show that X is 0.2-0.3x10"'° 
dyne cm sr"', in which Xm is about 0.15 - 0.2 x lO-"' dyne cm 
sr"' (Pizzo et al. 1983 1. Although the measured magnitude of 
X is smaller than that computed in the Weber-Davis model, the 
distribution of angular momentum flux between particles and 
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magnetic stresses is largely compatible with their prediction. 
An equally important finding concerns further distribution of 
X.P between two major ion species in the solar wind, namely, 
protons and alpha particles. Alpha particles are found to carry 
an angular momentum flux in the direction of counter-rotation 
with the Sun. This flux is substantial enough to offset the proton 
contribution which could be comparable to the magnetic one. 
This finding cannot be addressed by the Weber-Davis model 
where the solar wind was treated as a bulk flow and only pro- 
tons were considered. 

Apart from being essential in the problem of solar angu- 
lar momentum loss, the azimuthal ion motions may also pro- 
vide a possible means to resolve a long standing observational 
puzzle, namely that alpha particles are observed to flow faster 
than protons in the fast solar wind. The differential stream- 
ing in the fast wind could be as pronounced as 150 km s"^ 
at 0.3 AU before decreasing to about 40 km s"' at 1 AU 
(Marsch et al. 1982 1. Such a behavior has yet to be under- 
stood. One possible mechanism is the coupling between the az- 
imuthal and meridional motions, facilitated by the strong mag- 
netic field (McKenzie et al. 1979; Hollweg & Isenberg 1981 1. 
Although the Poynting flux may still be negligible (cf. Acuna 
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& Whang ^W!B, Alexander & de La Torre ( fT^ . Hu et 
al.(|2003 1), the proposed coupling is expected to limit, at least to 
a non-trivial extent, the ion differential streaming. As pointed 
out by HoUweg & Isenberg ( 1983 1, one shortcoming of the for- 
mulation of McKenzie et al. J1979> is that protons are privi- 
leged over alphas: The azimuthal magnetic field is assumed to 
be determined entirely by the protons whose azimuthal flow is 
neglected. Hence the formulation cannot properly account for 
the azimuthal dynamics of protons or alphas. In addition, the 
formulae are applicable only to the equatorial plane where the 
slow solar wind prevails at solar minimum. 

The goal of this paper is to extend the Weber-Davis model 
by including alpha particles self-consistently. This approach 
allows us to assess the individual contributions of ion flows 
and magnetic stresses to the angular momentum loss of the 
Sun. The effect of the coupling between azimuthal and merid- 
ional motions in limiting the proton-alpha diff'erential stream- 
ing will also be explored quantitatively. Given that the differ- 
ential streaming is more prominent in the fast wind, which in 
general flows out of the equatorial plane, it is necessary to for- 
mulate the model such that it treats both protons and alphas 
on an equal footing, and allows the solar wind to flow outside 
the equatorial plane. In this sense, this paper also extends the 
model of McKenzie et al. (« 1979> . 

The paper is organized as follows. The derivation of the 
governing equations is given in the appendix. Section|2ldetails 
the physical model and also describes the assumptions on the 
background poloidal magnetic field and the ion heating mech- 
anism. A general analysis is then given in section|3] Section^] 
presents the numerical results and the eff'ect of the solar rota- 
tion. In section|5] the main results are summarized. 



2. Model 

The solar wind is assumed to consist of electrons (e), pro- 
tons (p) and alpha particles {a). Since the role of alpha par- 
ticles is not necessarily minor, their contribution has to be self- 
consistently taken into account (Li et alj_997j. This is done by 
rewriting the momentum equations (Schunk [l977> in the flux 
tube frame, instead of the standard spherical coordinate system 
(r, 6, (p). Central to the derivation is that the ion-cyclotron fre- 
quency Q.k - ZkcBilmkC is many orders of magnitude higher 
than any other frequency present in the ion momentum equa- 
tions (McKenzie et al 1979) . Here is the mass of ion species 
k(k - p, a), Zk is the charge of species k in units of the electron 
charge e, Bi denotes the meridional magnetic field strength, and 
c is the speed of light. The derivation is provided in the ap- 
pendix, the resulting governing equations are quoted here to 
save space. 



2.1. Governing Equations 

As described in Appendix A, in addition to the impicit as- 
sumptions in deriving the five-moment transport equations 
(Schunk [19771 1. we make the following assumptions: 

1. Axial symmetry is assumed (d/d(f> = 0); 



2. The electron inertia is neglected in the electron momentum 
equation; 

3. Quasi-neutrality is assumed; 

4. Both viscosity and resistivity are neglected; 

5. Quasi-zero current is assumed, the only exception to this 
occurs when ion momentum equations are derived; 

6. The electric field in the magnetic induction law is con- 
vected by electrons, i.e., contributions like Hall eff'ects are 
neglected; 

7. The Spitzer law is used for the field-aligned electron heat 
flux, and the ion heat fluxes are neglected; 

8. We are only interested in steady state solutions. However, 
time-dependent equations are solved to yield steady state 
solutions. 

Given these assumptions, the governing equations take the 
form 
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where n^, Vs and denote the number density, velocity and 
temperature of species s (s = e, p, a), respectively. The species 
pressure is p^ - nJigTe,, where k^ is the Boltzmann constant. 
By assuming quasi-neutrality and quasi-zero current, we have 
Hg - Yik '^knk and Vg - Y^k ^k^k^klng {k - p, a). G is the gravi- 
tational constant, M© is the solar mass, and y = 5/3 is the adia- 
batic index. Coordinate / measures the arclength of the poloidal 
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magnetic field line from the footpoint at the coronal base. Both 
the heliocentric distance r and colatitude 6 are evaluated along 
the poloidal magnetic field. The cross-sectional area of the flux 
tube, a, scales as a oc O is the magnetic azimuthal an- 

gle, defined by tan<l) = B^/Bi. The poloidal magnetic field Bi 
and the heat deposition to ion species k, denoted by Qi^, will be 
specified in sections lT^ and E31 respectively. 

The energy and momentum exchange rates 6Es/6t and 
6 Ms /St are due to Coulomb collisions of species s with all the 
remaining ones (Schunk .l977j . 



to be apportioned between protons and alpha particles by 



-J— = 2_i"s'nsVst<l'st(vt - Vs), 



6t 

6t ^ nis + nil 



(8) 



[3kB(T^ - rs)^st + m^(v^ - Vs)'Ost] . (9) 



Expressions for the collision frequency Vst as well as correction 
factors Ost and ^'st have been given by, e.g., Li et al. fl997l l and 
will not be repeated here. In the computation, the Coulomb log- 
arithm In A is taken to be 2 L The electron thermal conductivity 
K, in Eq. Oil is 7.8 x 10"^ erg K''^'^ cm"' s"' (Spitzer[l962|l. 

2.2. Background poloidal magnetic field 

To avoid complications associated with the cross-field force 
balance, we choose to prescribe the background poloidal 
magnetic field by adopting an analytical model given in 
Banaszkiewicz et al. ( 1998"!. In the present implementation, 
the model magnetic field consists of dipole and current-sheet 
components only. A set of parameters M = 3.6222, Q - Q, 
K - 1.0534 and oi = 2.5 are chosen such that the last 
open magnetic field line is anchored at 6 - 40° on the Sun, 
and the poloidal magnetic field strength is 3. 3y at - 70° 
at 1 AU, compatible with Ulysses measurements (Smith & 
Balogh .1995 ). 

Figure shows the magnetic field configuration in the 
meridional plane. The thick solid line represents the field line 
along which we will find solar wind solutions. This field line is 
rooted at colatitude 3L5° on the Sun, and reaches 70° at 1 AU, 
which corresponds to the edge of the fast stream observed by 
Ulysses(McComas et al. 2000.) . Plotted in Fig.^i is the radial 
profile of the poloidal magnetic field strength Bi along the des- 
ignated field Une. 



2.3. Ion heating 

To produce fast solar wind solutions, an empirical energy flux, 
launched from the Sun and in the direction of B, is assumed to 
heat ions only. This energy flux is assumed to dissipate at a rate 
Q with a characteristic length i.e.. 



Bi I I 
Q - Fe exp 

BieIcI \ Id 



(10) 



where Fe is the input empirical flux scaled to the orbit of the 
Earth, Re - 215 Rq, Rq being the solar radius. Moreover, Bie 
is the poloidal magnetic field strength at Re- Q is then assumed 
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1 + A 



1 



1 +A 



Pp 



(11) 



where pi - ni^m^ (k - p, a) is the ion mass density, and x is 
a parameter indicating the degree by which the alpha particles 
are preferentially heated, with ^ = 1 standing for the neutral 
heating: heating rate per ion is proportional to its mass. 
In the computations, the following parameters 



Fe^I.S erg cm-2 s-\ Id = 1.35 Rq, 

Xc = 1.5 



X = — tanh 



0.3Rz 

are chosen to yield a fast solar wind solution. As can be seen, 
X varies smoothly from Xc in the inner corona to 0.8 far from 
the Sun with a rather steep transition occurring at 5 Rq. A pref- 
erential heating that favors alpha particles in the inner corona 
iXc > 1) is necessary to produce a positive relative speed 
Vol -Vpi- 



3. General Analysis 

Before solving equations Q to to find solar wind solutions, 
one can conduct an analysis to reach some general conclusions. 



3.1. Alignment conditions 

Equation ^ derives from the (p component of the magnetic 
induction law. For a steady state, it can be integrated to yield 



Bd, 

- Orsmfl = —Vei- 
Bi 



(12) 



The constant of integration Q. can be identified as the angular 
rotation rate of the flux tube. Combining Eq. J12l i and Q , one 
finds 



Bit, 

Vp0 - Qrsmfl = — V, 
Bi 



B(j, 

Va4, - Or sine - —v„i. 

Bi 



(13) 



That is, in the frame strictly corotating with the Sun, all species 
(electrons, protons and alpha particles) flow along the mag- 
netic field. The alignment conditions were first recognized by 
Parker ( il958< . and have been extended to general axisymmet- 



rical MHD flows by, e.g.. Low & Tsinganos J 19861 1 and Hu et 
al. ll2003l . 

3.2. Angular momentum conservation law 

In a steady state, equation ^ leads to 
BaBi 



rsinO 



^>0 



Pa^al 

+ Vo 

PpVpl 



4-7TppVpl 



(14) 



where the tube invariant L comes from the integration. The 
physical meaning of L can be better seen by noting that the 
constant 



-C - ppVpiL — R^ - £p + Xtr + -Cm 



(15) 
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is the angular momentum loss per solid angle, where 



3.4. Energy conservation law 



Zk - r&inOpkVkiVk^ — R\, -Cm = sin 



^B^Bi a 

47T Ue' 



Ri 



(16) 



with k - p,a. Subscript E denotes values evaluated at 
Re - I AU. Obviously, both outflowing particles and magnetic 
stresses contribute to the angular momentum flux. 

The conservation law for the angular momentum, Eq. lll4> . 
is valid for an arbitrary flux tube in an azimuthally symmetric 
solar wind. The single-fluid version (or equivalently the two- 
fluid one) of this conservation law has already been obtained 
by, e.g.. Low & Tsinganos (1986^ and Hu et al. (l2003t . 

3.3. Expressions forvp^, Vq,^ and 

The alignment condition, Eq. ( I13t . together with Eq. (I14t leads 
to 
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Qr sin 
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where Mj, Mp and M^, are defined as 

,.2 



Mj=Ml+ Ml, 
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BjlAnpk 



(17) 



(18) 



(19) 



(20) 



with k - p,a. The poloidal Alfvenic Mach number Mr is thus 
comprised of both Mp and Ma- 

For the solar wind, Mj 1 is valid near 1 Rq, but Mr » 1 
holds at 1 AU. Hence, there must exist a point between 1 Rq and 
1 AU where Mj = 1 . At this location, which will be termed the 
Alfvenic point, is singular unless the numerator in Eq. jl9t 
vanishes. 



1 + 



PaVaO 
PpVplI' 



Qri sin^ e„ 



(21) 



where subscript a denotes values at the Alfvenic point. We have 
employed the fact that the ion mass flux ratio PaVai/PpVpi is 
a constant. The angular momentum loss per solid angle then 
becomes 

£ = MQrl sin^ da, (22) 
where 

M - (ppVpl + paVal) — Re 

is the mass loss rate per solid angle of the solar wind. 

Hence the conclusion of Weber & Davis J 19671 1 still holds: 
the magnetic field helps the coronal plasma to achieve an ef- 
fective corotation to the Alfvenic point, as long as the poloidal 
Alfvenic Mach number Mj is defined by Eq. ( I20> . 



Combining the governing equations in the steady state, one can 
derive an energy conservation law. 



a 

Qe 



PpVpl- 



+ PaVal- 



(ppVpl + PaVal) —GMq ( - - I 

^ ' Qe \Rq rj 

a B^ / X 
— —\v„B^-v,^Bi) 

Y [peVel + PpVpl + PaVal) 



a B. 

Ue 4-7T 

a y 



ClE y 

-^.rf^^cos^O 
fl£ al 



r 



— {Qp + Qa)Al ■ 

Ue 



r. 



(23) 



where the constant T is the total energy flux scaled to Re- The 
terms on the left hand side of Eq. ( I23> correspond, respectively, 
to the kinetic and potential energy fluxes, the Poynting flux, the 
enthalpy flux, the electron conductive flux, and the source term 
due to the heat deposition. The ratio of the Poynting flux to T 
will be used to assess the relative importance of the Poynting 
flux in the energy budget. 

4. Numerical Results 

Equations ([0 to are solved by using a fully implicit numeri- 
cal scheme (Hu et al. 1997 1. From an arbitrary initial guess, the 
equations are advanced in time until a steady state is achieved. 
The computational domain extends from 1 Rq to 1.2 AU. At 
1 Rq, ion densities as well as species temperatures are fixed. 



rip = 1.5 X lO'^cm" 



Te = Tp =Ta^ lO^K, 



(ua/np) = 0.06, 



while Vpi and Vai are specified to ensure mass conservation. Vp^ 
and Bif, are evaluated in accordance with equations (I17> and 
(I19> . where L is computed at the grid point immediately adja- 
cent to the base. At the outer boundary (1.2 AU), all dependent 
variables are linearly extrapolated for simplicity. We also take 
Q - 2.865 X 10"^ rad s"'. For the steady state solutions pre- 
sented in this paper, the maximum relative errors in the con- 
served quantities are smaller than 1%. 

Figure 12 displays the radial distribution of (a) the species 
densities n^, rip and «„, (b) poloidal flow speeds Vpi and Vai, 
and (c) species temperatures Tp and Ta- The model yields 
the following parameters at 1 AU, 



ripVpi - 2.3 X lO^cm ^s Vpi - 660 km s 
v„, - Vp/ = 49 km s"' , ria/rip = 0.0445 



which agree well with the Ulysses observations of the fast 



wind (McComas et al. l2000l l. In addition, the modeled electron 
density fits observations reasonably well in the inner corona 
(Fig.|2t). However, without considering the non-thermal con- 
tribution, the modeled proton temperature is higher than that in- 
ferred from UVCS measurements (Fig.|2l;). Moreover, Tp and 
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Ta at 1 AU are much smaller than values given by in situ mea- 
surements (McComas et al. I2000> . The poor match is due to 
the oversimplified heating function. As for the poloidal flow 
speeds, the alphas initially fall behind the protons below 2 Rq 
beyond which a positive Av/ - Vai - Vpi develops. Vai reaches a 
maximum around 66 Rq, and starts to decrease thereafter 

To examine the differential streaming further, Av; = Vai-Vpi 
is plotted in Figure|3t. The poloidal flow speeds of protons ( vpi) 
and alpha particles (vv/) are replotted in Fig. O) (a different 
scale is used, see Fig.|2j5). In addition, model results from the 
corresponding computation that neglects the solar rotation (i.e., 
Q s 0) are plotted as dotted lines for comparison. For the ease 
of description, we shall call the model with (without) azimuthal 
components the 1 .5-D ( 1 -D) model. It is found that the effect of 
the azimuthal components on the poloidal dynamics can be ad- 
equately represented by the flow speed profiles. Below the local 
maximum of 78.6 km s"' at 7.28 Rq, Fig.|3t shows no differ- 
ence in the Av; profile between 1-D and 1.5-D models. The dif- 
ferential streaming, Av;, for both models plummets from nearly 
zero at the coronal base to about -44.6 km s"' at 1.44 Rq, and 
rises thereafter to the local maximum. Interestingly, in the 1-D 
model, beyond the local maximum Av/ undergoes only a mod- 
est decrease to 66.3 km s ' at 1 AU, while in the 1.5-D model 
Av/ is 48.7 km s"' at 1 AU. This further reduction in the dif- 
ferential streaming is achieved through a slight rise in the Vpi 
profile accompanied by a modest deceleration of alpha parti- 
cles (Fig.H)). 

This behavior is not surprising since in the poloidal mo- 
mentum equation (Eq. (|2}), 



VkV 



dvki 
dl 



2 ^ , . 



. dvk4, d , 

tan Ovt) h vij, — In r sm ( 

' dl ^dl 



d_ 
di 



^sec2. 
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d qV- 



^ sin^ ( 



dl 



(24) 



can be obtained when the alignment condition, Eq. ( I13t . is 
used. When viewed in the frame corotating with the Sun 
(vjt/ sec O is the ion speed seen in that frame), the solar rotation 
ensures that all particles move in the same centrifugal potential 
(Q^r^ sin^ QjT). Neglecting all other contributions, and taking 
the difference of the proton and alpha version of Eq. (I24l i. one 
arrives at 



dl 
or 

ivli - vl,) oc cos2 (D. 



(25) 



With the development of the magnetic azimuthal angle, cos^ O 
decreases monotonically with increasing distance (cf. Fig.|^). 
As a consequence, the differential streaming Av; decreases. 
Figure|3]can be seen as a direct illustration of the effect of solar 
rotation in limiting the ion differential streaming, predicted by 
McKenzie et al.l fDT^ and Hollweg & Isenberg ( flWTt . 

Figure 0] displays the radial profiles of (a) - tan O = 
-B^/Bi, (b) the azimuthal speeds of protons (v^^), alpha par- 
ticles (Va^) and electrons (v^^), and (c) the specific contribution 



of protons (^p - JLpjX,), alpha particles {^^ - -Cal -O and the 
magnetic stresses {^m - JLmI-O to the angular momentum flux 
(cf. Eq. M5V). In addition, the sum ^p + f^,, which gives the 
overall particle contribution ^p, is also plotted. Given in dotted 
line is the ratio of the Poynting flux to the total energy flux 
(cf. Eq. (I23» . In Fig.|3;, the dashed line is used to plot nega- 
tive values. The asterisks in Fig.|4j) denote the Alfvenic point, 
which is located at r,, = 11 .8/?©. 

From Fig. |4t, it is obvious that only beyond, say 10 ^o, 
does a spiral angle O develop. This can be explained in view of 
equations (I12> . J13> : within 10 the left hand side is much 
smaller than the poloidal flow speed on the right hand side 
for any species. On the other hand, in interplanetary space, 
the species azimuthal speed is much smaller than Qr sin 9, the 
Parker theory for the spiral magnetic field is recovered, i.e., 
tanO - B^/Bi x -Q.rsm6/vi, where v/ can be taken as the 
poloidal speed of any species. 

In the inner corona, both protons and alpha particles tend 
to corotate with the Sun: Vp^ and v„0 are positive (Fig.|?J)). The 
azimuthal speed of the alpha particles is slightly larger than that 
of the protons below 2 Rq, and from there on, the alpha particles 
are gradually turned opposite to the solar rotation. becomes 
negative beyond 5.71 Rq, eventually reaches -24.7 km s"' 
at 1 AU. On the other hand, the proton azimuthal speed Vp^ in- 
creases from a local minimum of 0.78 km s"' at 12.5 Rq mono- 
tonically to 4.8 km s"' at 1 AU. 

The behavior of the azimuthal flow speeds can be explained 
by Eqs. (tTTt and (tTsl . Near the coronal base, both Mj, and 
are far from unity. It then follows from Eqs. (I17> and (I18> that 



Qr sin 6, vv^ a; Qr sin 0. 



(26) 



At r < 2 Rq, a negative Av; makes Vaip slightly larger than 
Vp^. For r > rg, the solar wind expands almost radially. As 
a result, Mj/va - {AnpkVkilBifjBi cc (k - p,a) holds 
fairly accurately. The variation of Vpi beyond rg is very mod- 
est. We therefore have » 1 for r » Tq. From the identity 
^a/^p - ^"^'ai/vpi, it follows that is a substantial fraction 
of Mjj given that the ion mass flux ratio 77 = PaVai/PpVpi is 0.19 
in this solution. Hence close to 1 AU, the azimuthal speeds of 
both protons and alphas are determined by the terms associated 
with the differential streaming in Eqs. (I17t and (I18> . namely, 
for r » Tfl, 



Or sin 6 



Val - Vpi 



-Qr sin 6 



1 + rjVallVpi Vpi 

1 Val - Vpi 

1 + rjVailVpi Vpi 



(27) 



As a result, Vp^lva,p ~ -t] holds. However, this asymptotic be- 
havior of the ion azimuthal speeds for r ^ rg does not hold in 
general. If the alpha abundance is far from unity, the azimuthal 
magnetic field will be solely determined by protons, and Vp^ 
should behave like Vp^ oc r ' for r :s> when the differential 
streaming term in Eq. is neglected. 

Now let us move on to Fig.|3;. It can be seen that, from the 
coronal base to 1 AU, magnetic stresses play a dominant role 
in the total angular momentum budget, the particle contribu- 
tion is no more than 2.6%. However, the individual angular 
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momentum flux carried by protons or alpha particles is not nec- 
essarily small in magnitude. As a matter of fact, protons con- 
tribute more to the total angular momentum flux than magnetic 
stresses do beyond 101 Rq- However the proton contribution 
is virtually canceled by the alpha particles that counter-rotate 
with the Sun. This can be understood in light of equations ( I17> 
and ( I18> . As has been described, far away from the Alfvenic 
point, r » Tfl, both Vp^ and are mainly determined by the 
terms associated with the diff'erential streaming. From the iden- 
tity PaVaiMjjvpi = PpVpiM^/vai, One can see that for r » r,,, ^„ 
and ^p tend to have the same magnitude but opposite sign. At 
this point, we can also see from the dotted line in Fig.^ that 
although the solar rotation introduces appreciable difference 
in the meridional dynamics, the Poynting flux never exceeds 
0.12% of the total energy budget. Needless to say, its contribu- 
tion to the solar wind acceleration is in fact determined by its 
difiference between 1 Rq and 1 AU. 

At 1 AU, the model yields a total angular momentum loss 
oi £. = 0.17 X lO^*' dyne cm sr in which the magnetic part is 
Xm = 0.165x10^" dyne cm sr' ' , consistent with measurements 
(Pizzo et aL "19831 Marsch & Richter [19841. However, the ab- 
solute azimuthal speed Vp^ or Vq,^ is larger than the measured 
values (although these quantities can only be determined with 
a modest precision). Moreover, never turns negative, in this 
sense at variance with the measurements: Pizzo et al. ( 1983]) 
and Marsch & Richter ( 1984i showed that particles in the fast 
wind tend to carry a negative angular momentum flux. Pizzo et 
al. 11983 1 suggested that the discrepancies between the model 
and measurements may be removed by including the stream in- 
teraction in the super- Alfvenic region. This is however beyond 
the scope of this paper 

5. Concluding remarks 

The main aim of this paper is to extend the Weber-Davis anal- 
ysis (Weber & Davis 19671 1 on the transport of the angular 
momentum from the Sun by including alpha particles and by 
allowing the solar wind to flow out of the equatorial plane 
in an axisymmetrical configuration. Following McKenzie et 
al. < 19791 1. we exploit the fact that the gyro-frequency of ions is 
many orders of magnitude higher than any other frequency in 
ion momentum equations. From this it follows that the differ- 
ence between proton and alpha velocities must be in the direc- 
tion of the magnetic field. Using this alignment condition, the 
governing equations are then derived from the standard five- 
moment transport equations. 

The model equations also enable us to examine quanti- 
tatively the effect of azimuthal components in limiting the 
proton-alpha differential streaming in the fast wind. For sim- 
plicity, we choose to solve the governing equations on a 
prescribed poloidal magnetic field line located at a colati- 
tude of 70° at 1 AU, corresponding to the edge of the fast 
stream observed by Ulysses at solar minimum conditions 
(McComas 2000). The effects of the azimuthal components 
on the meridional dynamics, if any, are optimal in this regard. 
These effects are directly shown by a comparison of two mod- 
els with and without azimuthal components. 

The main results can be summarized as follows: 



1 . The general analysis concludes that, in agreement with the 
Weber-Davis model, the magnetic field helps the coronal 
plasma to achieve an effective corotation from the coronal 
base to the Alfvenic radius, where the poloidal Alfvenic 
Mach number Mj - 1 . Mr has to include the contribution 
from alpha particles (Eq. j20ll). 

2. In the low latitude fast solar wind, the angular momen- 
tum loss from the Sun is almost entirely due to magnetic 
stresses. The proton contribution, which can be as impor- 
tant as the magnetic one in interplanetary space, is offset 
by alpha particles that develop an azimuthal speed in the 
direction of counter-rotation with the Sun. 

3. The Poynting flux associated with the azimuthal compo- 
nents is negligible. Nevertheless, the solar rotation has an 
appreciable effect in limiting the proton-alpha differential 
streaming in fast solar wind streams at low latitudes in in- 
terplanetary space. 

Although the fast solar wind solution is largely compatible 
with in situ measurements in terms of the ion mass fluxes and 
terminal speeds, it fails in a detailed fashion. For instance, the 
model is not able to predict a proton temperature profile con- 
sistent with UVCS measurements in the inner corona, nor does 
it predict an ion differential speed as large as 150 km s"' at 
0.3 AU to be comparable with the Helios observation (Marsch 
et al.T5ST). Hence, including the azimuthal components can- 
not solely account for the deceleration of alphas relative to 
protons in interplanetary space. More sophisticated mecha- 
nisms, the ion-cyclotron resonance for instance, are expected 
to alleviate the discrepancies (e.g., Li .2003 1, but can hardly 
help achieve a satisfactory match (e.g., Hu & Habbal ll999> . 
Nevertheless, such a direction is for sure worth pursuing and is 
left for a future study. 

The model also suffers from the inconsistency that the force 
balance in the direction perpendicular to the poloidal magnetic 
field is replaced by prescribing a background magnetic field. In 
a more rigorous treatment, the poloidal magnetic field should 
be derived self-consistently. In principle, such a task can be 
accomplished by adopting an iterative approach: the parallel 
and perpendicular force balance are solved alternately until a 
convergence is met (Pneuman & Kopp [T97Tl Sakurai ll985> . 
By doing so, the angular momentum loss from the Sun can be 
obtained self-consistently for all poloidal flux tubes. An accu- 
rate estimate of the duration over which the angular momentum 
of the Sun is completely removed is then possible (see Hu et 
al. l2003l . 

The present paper is aimed at presenting a rather general 
analysis of the angular momentum loss from a magnetized ro- 
tating object for flows assuming axial symmetry and incorpo- 
rating two major ion species. Although for the present Sun, the 
centrifugal and magnetic forces are so weak that they have lit- 
tle impact on the meridional dynamics (especially below the 
Alfvenic point), a similar study as presented in the text can be 
carried out for stars that rotate at a faster rate or have a stronger 
magnetic field than the Sun. 
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Fig. 1. (a): The poloidal magnetic field configuration given as contours of the magnetic flux function. The equator points upward. 
The line of force on which the model equations are solved is displayed by the thick contour This field line is located atO - 70° at 
1 AU and originates from 3L5° on the Sun. (b): Radial distribution of the poloidal magnetic field strength Bi along the designated 
field line. At 1 AU, Bi is 3.3y. 




1 10 100 1 10 100 1 10 100 

^Ro) KRo) ^(Ro) 



Fig. 2. Results derived from a 1 .5-D 3-fluid solar wind model which incorporates the azimuthal components self-consistently. The 
radial distribution of (a) the densities of protons Hp and alpha particles n„ (solid lines), as well as electrons tig (dashed line), (b) the 
poloidal flow speeds of protons (vpi) and alphas (vai), and (c) the temperatures of electrons (Te), protons (Tp) and alpha particles 
(Ta). The error bars in (a) are the upper and lower limits for the electron density derived by Fisher & Guhathakurta ( 1995 1. The 
error bars in (c) represent the uncertainties of UVCS measurements for the effective proton temperature reported by Kohl et 
al. (I1998> . Please note that both measurements are made for polar coronal holes. 
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Fig. 4. Radial distributions of (a) - tan (5 = -B^jBi where C) is the magnetic azimuthal angle, (b) the azimuthal speeds of 
protons Vp^, alpha particles v„0 as well as electrons v^.^, (c) the relative importance of the proton fluid - Upl H, the alpha fluid 
- -Lai -L, the sum of the two - ^p + ^a, and the magnetic stresses - -Cm/-C in the total angular momentum loss of the 
Sun (please see Eq. (I16> in text). In addition, the ratio of the Poynting flux to the total energy flux, ^, is plotted as dotted line. The 
dashed line represents negative values. In panel (b), the asterisks denote the Alfvenic point, where the poloidal Alfvenic Mach 
number (defined by Eq. (I20» equals unity. 
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Appendix A: Derivation of the governing 
equations 

In this appendix, it is shown how the 5-moment transport equa- 
tions are reduced to the governing equations in section lzTI The 
approach adopted here closely follows that by McKenzie et 
al. (I1979> (see also Hollweg & Isenberg |l981> . The original 
derivation of McKenzie et al. f 1979^ is restricted to the equa- 
torial flow, and ions other than protons are treated as test par- 
ticles. Employing the same spirit, we extend their derivation 
to general flows assuming axial symmetry. In addition, all ion 
species are treated on an equal footing, which is particularly 
important for the solar wind since alpha particles can not be 
seen as test particles. The central point is that, due to the pres- 
ence of a strong magnetic field (in the sense that the ion gyro- 
frequency is many orders of magnitude higher than any other 
frequency in the momentum equations), the diff'erence vector 
between proton and alpha velocities must be ahgned with the 
magnetic field. 

A.1. General momentum equation 

First of all, let us examine the momentum equation for species 
s (Schunk[T973, 



nsifis 



dvs 
'dt 



+ Vs • VVc 



+nsms — Tj—r - n^e^ 



E+ -v^xB 



6t 



= 0. 



(A.l) 



As usual, species s is characterized by its density n^, velocity 
Vs, mass nis, electric charge es and pressure ps- e^ can also 
be measured in units of electron charge e, i.e., es = Zse with 
Ze = -1 by definition. The momentum exchange rate 6Ms/6t 
is due to the Coulomb frictions. It is customary to neglect the 
electron inertia (nig = 0). As a result, the electrostatic field E 
can be expressed as 



E = --Ve 

c 



xB 



_!_ + 



dt 



(A.l) 



(A.3) dB^ 
dt 



Substituting the expression for E into the magnetic induction 
law 

— H- cV X £ = 0, 

dt 

one then arrives at 
dB 

— - V X (v, X B) = 0, 
ot 

where B is the magnetic field. The terms in Eq. ( IA.2> other than 
the motional electric field -Ve x B/c are many orders of mag- 
nitude smaller and thus have been neglected. All terms have to 
be kept when E is substituted into the ion momentum equation 
however. 

The plasma in question consists of two ion species, protons 
(p) and an additional one (/). (Subscript / is used here to indi- 
cate that in principle the equations to be developed are also ap- 
plicable if other ion species than alpha particles is considered.) 



As the frequency in question is well below the electron plasma 
frequency, the expression for n^ follows from quasi-neutrality, 

Hp = n,, + Z/H/. (A. 4) 

Neglecting the displacement current in the Ampere's law, one 
finds the expression for Ve 

_ npVp + ZinjVi ^ j 
«(, n^e 



(A.5) 



where j - cV x B/4jt is the electric current density. 

Substitution of Eq. (IA.2> into Eq. iA.l\ for ion species k 
(k - p, i) then leads to 

dVk Vpk Zklpe 

— +Vk-lVk + + 

at nicnik ngnik 



GMs 
1 

nkntk 



Annenik 



(VxB)xB 



5Mk Z^nt SMe 
+ - 



St 



St 



,^!!i£^(v,-v.)xB.O, 

nikc ne ^ ' 



(A.6) 



where subscript j stands for ion species other than k, namely, 
i - p for k - i and vice versa. Note that, when deriving equa- 
tion ( IA.6> . we have used equation (IA.5t to evaluate the electron 
velocity in the expression for E. The electric current j can be 
dropped when Ve is evaluated elsewhere. This is because, in the 
context of the solar wind, j is the large-scale electric current 
and is negligible since the spatial scale at which the magnetic 
field evolves is well beyond the proton inertial length. 

A.2. Alignment conditions for electrons 

Now it becomes necessary to separate explicitly the poloidal 
and azimuthal components of the magnetic field and species 
velocities, namely, 

B - Bp + B^^, Vs = Vsp + Vs00, 

where subscript P stands for the poloidal component. The as- 
sumption of azimuthal symmetry (d/d<p-Q) allows Bp to be ex- 
pressed in terms of the magnetic flux function t//(r, 0; f), i.e.. 



B, 



Vi// X 



rsinO 

The magnetic induction law, Eq. iA.3\ . can then be rewritten as 



dip 

^ + v,p ■ Via = 0, 
ot 



+ r sin 0V ■ 



1 n W V ' 



r sinfl 



= 0. 



(A.7) 
(A.8) 



For a steady state, from Eq. ( IA.7> follows Vep ■ V^/r - 0, which 
is equivalent to 



VeP XBp -Q. 



(A.9) 



In other words, the poloidal components of the electron ve- 
locity and magnetic field are strictly parallel. In light of this 
alignment condition, equation (IA.8> can be shown to reduce to 
Eq. (|6} in the text. 
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A.3. Equations cast in tlie flux tube frame 

It proves useful to work in the flux tube frame, whose base 
vectors are 

ei — Bp/Bp, e^ — 4>,^2 — e^ Xei. 

By definition, the magnetic field has only two components, i.e., 

B - Biii + £363, 

We restrict ourselves to time-independent axisymmetrical 
flows only. Taking the dot product of e^ with Eq. ( IA.6> . one 
arrives at 

VkidiVki - vl^di Inrsinfl 

+Vj2ei ■ (e2 ■ Ve2) + Vk2d2Vki + vuVkiei ■ (ei ■ Ve2) 

1 Z, GMo 

H Oipk H Oipe H 0\ In r 

nknik n^nik r 

2k 

+ — B3 {dxB^+B^di In r sin 6*) 

A-menik 

nj B3 Zjnj 

coivji - Vki) + ^ki- (vj2 - Vk2) = 0. (A.IO) 

Akne Bi n^ 

Similarly, taking the dot product of Eq. iA.6\ with 62 and 63 
results in, respectively. 



^ - v^3<92lnrsin( 



+Vk2d2Vk2 + Vk\d\Vk2 + nyVkih ■ ih ■ Vei) 

1 Z, GMo 

+ 02Pk + OlPe + 02inr 

n^iik n^nik r 



Zk 



AnneMk 
n 



do 



B^,+Bl 



+ B\d2 In r sin 6 



--r—CQ{vj2 - Vk2) 

Akn, 



(V;3 - Vk-i) - (Vjl - Vk\) 



(A.ll) 



and 



Vk\ (dyVki + Vk-idx Inrsinfl) 

+Vk2 (d2Vk3 + Vk3d2 In r sin 0) 
Zk 

-- — —Bi (diBi + B3di In r sine) 
Annemk 

rij 
Akn, 
Zjnj 

-^kl - Vk2) = 0, 



(A.12) 



ZkcBilnikC is the gyro-frequency for species k. d„ - e„ ■ V is 
the directional derivative operator along e„ (n = 1, 2, 3). 

The ion gyro-frequency is many orders of magnitude higher 
than any other frequency in the momentum equation. This has 
two consequences. First, from equation (IA.12> . Vj2 - Vk2 is far 
smaller than Vk-i from an order-of-magnitude estimate. Since 
Ve2 - (see Eq. (IA.9» . both Vp2 and va should be very small 
and can be safely neglected unless they appear alongside the 
ion-cyclotron frequency. Second, Eq. JA.IU leads to 



V,-3 



B3 



^(v/i-vpi). 



(A.13) 



In other words, the ion velocity difference is aligned with the 
magnetic field. This is Eq. 0. 
Solving Eq. 



( IA.12> for Vj2 - Vk2 and then substituting it 
into Eq. ( lA.lOL one arrives at the poloidal momentum equa- 
tion (Eq. It is interesting to note that the magnetic field 
does not appear explicitly in this equation (except for the term 
tan ft), although it plays an essential role in coupling the az- 
imuthal and meridional motions. Combining Eq. ( IA.12> for p 
and /, one can obtain Eq. (|5}- The partial differentiation with re- 
spect to time t in the equation is merely for numerical purpose. 

In closing, we note that the p and / versions of Eq. (lA.lH 
can be combined to yield a force balance condition in the 62 
direction. 



- Ppv],T,d2 In r sine 
PiV%d2 In r sine 



+ d2Pp + d2Pi + d2Pe + (Pp+ Pi) 92 In r 



1 

4^ 



d2 



b] + bI 



+ B\d2 In r sin 



= 0. 



(A. 14) 



If further expressing the geometrical coefficient the differ- 
entiation 82 and the Lorentz force in terms of the magnetic 
flux function i/r, one can eventually derive a second-order quasi- 
linear partial differential equation (PDE) for i/r. This PDE can 
then be solved by using the approaches proposed by Pneuman 



& Kopp ( IT97Tb or Sakurai ( [19851 



where 

l/e2-(ei ■ Vei) 

is the (signed) curvature radius of the poloidal magnetic field 
line, while 



ne 



Co = — Vpe<S>pe + —Vpi^pi + V/gOie 

n,- ne 



Zfri! 
n„ 



is a coefficient associated with Coulomb frictions. Here A,- = 
nii/nip is the mass number of species In addition, Q.ki = 



